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摘 要 :针对 地 下 目标 体 周围 介质 参数 估计 问题 ,提出 一 种 基于 图 像 炉 的 介质 参数 估计 方法 。 首 先 推导 了 介质 条 件 下 频率 

分 集 阵列 探 地 雷达 的 点 扩展 函数 ,分 析 了 点 扩展 函数 与 后 向 投影 成 像 算法 之 间 的 联系 ,然后 将 频率 分 集 探 地 雷达 与 超 宽带 

探 地 雷达 成 像 结 果 做 了 对 比分 析 。 对 于 同一 块 包含 目标 的 待 成 像 区 域 而 言 ,介质 参数 对 电磁 波 的 传播 速度 有 较 大 的 影响 ， 

当 采 用 不 同 的 介质 参数 分 别 对 该 区 域 进行 成 像 时 ,得 到 的 成 像 结果 亦 不 相同 。 对 不 同 参 数 下 的 成 像 结 果 进 行 图 像 依 的 计 

算 , 成 像 结果 的 图 像 粹 值 越 小 , 则 在 使 用 该 参数 时 成 像 的 聚焦 度 越 好 ,所 对 应 的 介质 参数 越 接 近 其 真实 值 。 实 验 结 果 表 明 ， 
在 介质 的 电导 率 不 为 零 时 ,频率 分 集 探 地 雷达 较 超 宽带 探 地 雷达 在 目标 定位 成 像 方面 具有 更 好 的 效果 , 且 当 目标 为 细 长 目 
\ 体 时 ,所 提出 的 参数 估计 方法 能 够 对 目标 体 周 围 介质 参数 进行 有 效 估计 。 
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stract: Aiming at the problem of medium parameter estimation around buried objects, a method of medium parameter es- 


astimation based on image entropy is proposed. Firstly, the point spread function (PSF) of frequency diversity array ground 


etrating radar (FDA-GPR) in medium is derived, and the relationship between PSF and back projection imaging algo- 
置 hm is analyzed. Then, the imaging results of FDA-GPR and UWB-GPR are compared and analyzed. For the same region 
containing the target to be imaged, the medium parameters have a great influence on the propagation velocity of electromag- 
netic wave. When different medium parameters are used to image the region, the imaging results are also different. The im- 
age entropy of the imaging results under different parameters is calculated. The smaller the image entropy of the imaging 
results is, the better the focusing degree of the imaging is, and the closer the corresponding medium parameters are to their 
true values. The experimental results show that: when the conductivity of the medium is not zero, FDA-GPR has better per- 
formance in target location and imaging than UWB-GPR, and when the target is a slender target, the proposed parameter 
estimation method can effectively estimate the parameters of the medium around the target. 
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在 工程 物探 领域 , 探 地 雷达 (ground penetrating 传统 的 探 地 雷达 大 多 采用 的 是 超 宽带 信号 ,即时 域 窗 
radar, 简 称 GPR ) 是 使 用 最 为 广泛 的 方法 之 一 。 脉冲 信号 。 然 而 窄 脉冲 信和 号 的 瞬时 带宽 过 大 ,不 仅 会 
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给 发 射 机 前 端的 硬件 实现 带 来 困难 ,使 成 本 过 大 ,而 
且 地 下 介质 也 可 能 会 存在 损耗 ,使 得 超 宽带 信号 产生 
畸变 , 易 造 成 目标 物理 特征 解 译 不 精确 。 对 此 ,可 
采用 多 个 鹤 带 信号 合成 宽带 的 形式 ,如 频率 分 集 阵 列 
(frequency diverse array, 简称 FEFDA)5 ,以 避免 
GPR 发 射 超 宽带 信号 时 带 来 的 大 带宽 和 非 线性 的 问 
题 , 从 而 降低 数据 的 宛 余 度 ,提高 信号 处 理 的 实时 性 
与 可 靠 性 

雷达 成 像 技 术 可 以 对 目标 的 位 置 和 轮廓 有 一 个 
大 体 的 呈现 , 现 有 成 像 算 法 主要 有 衍射 层 析 (diffrac- 
tion tomography, 简 称 DT)、 时 间 反 转 (reverse time 
migration, 简称 RTM) 和 后 向 投影 (back projection， 
简称 BP) 等 算法 中 。 其 中 BP 算法 可 通过 计算 阵 元 
到 成 像 区 域 像 素 点 的 距离 ， 从 而 对 电磁 波 在 介质 分 界 
面 处 发 射 的 折射 进行 精确 补偿 ,具有 很 大 的 实用 
性 倍 ). 在 对 目标 进行 成 像 时 ,需要 介质 的 相关 参数 ， 

有 的 参数 测量 通常 都 是 通过 抽取 部 分 样本 来 估计 
整 不 探测 区 域 背景 的 介 电 属性 ,这 种 方法 不 但 对 探测 
和 畏 员 有 有 破 环 性 ,而 且 效 率 较 低 。 为 此 ,引入 图 像 处 
理 领 域 中 图 像 箭 的 概念 ,通过 寻求 最 小 图 像 糯 值 ， 

和 进而 反 扒 出 介质 的 相关 电 特 性 参数 。 


硕 率 分 集 探 地 雷达 信号 模型 


在 远 场 条 件 下 ,频率 分 集 阵列 可 通过 波束 方向 图 
进 和 表征 。 而 对 于 近 场 的 情况 ,通常 采用 点 扩展 
消 数 (point spread function ,简称 PSF) 进 行 分 析 , 点 
扩 健 函数 可 表示 成 像 系 统 对 点 目标 的 响应 ,模型 示意 
图 项 图 1 所 示 。 


(Xp 5p) 
p» “Fr 


图 1 频率 分 集 阵列 探 地 雷达 模型 示意 图 


考虑 N 个 收发 同 置 阵 元 , 阵 元 间距 为 d ,各 阵 元 
发 射 窄带 信号 ,第 ”个 阵 元 发 射 的 信号 为 
S, (1) =A(z)exp(j2rf,t), (1) 
其 中 ,A(z) 为 包 络 , 罕 带 情况 下 ,其 变化 是 缓慢 的 。 
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言 号 频率 可 表示 为 
fi =fotmAfyrn=0ly "Ne—1, (2) 

其 中 :Af 为 单位 频率 增 量 ; f。 为 初始 载 频 ;mm 为 
NX1 频率 增 量 系数 向 量 的 第 n 个 元 素 ,频率 增 量 系 
数 向 量 m 由 0~N 一 1 的 计 个 整数 随机 排列 得 到 。 

假设 一 点 目标 P 位 于 介质 二 中 ,其 位 置 坐标 可 


表示 为 (zx,,z,)。 电 磁 波 在 介质 一 中 的 传播 速度 为 
cl 在 介质 二 中 的 传播 速度 为 c, ,第 n 号 阵 元 至 折射 
点 距离 为 7 ,折射 点 至 目标 点 距离 为 7,,,。 根 据 折 
射 定律 
sin 01 可 
sin0， cy” 
有 
(CI 三 人 
(Xp — Tr ca “0 
其 中 ,rj, .rs, 可 由 勾 股 定理 
ri 一 VCD 一 三 十 (一 过 
ran =V (xy xX,)” CR 之 


得 到 。 

联 立 式 (4) 和 式 (5) ,可 得 到 一 个 关于 折射 点 横 坐 
标 zx 的 四 次 方程 ,对 该 方程 进行 求解 ,可 得 到 4 个 
根 。 其 中 2 个 根 为 复数 , 因 折 射 点 的 坐标 必然 为 一 实 
数 , 故 对 这 2 个 根 不 予 考虑 , 男 外 2 个 根 均 为 实数 ,此 
时 观察 折射 点 的 位 置 ,发 现 该 点 处 于 发 射 阵 元 的 横 坐 
标 与 目标 点 的 横 坐 标 之 间 , 故 可 取 对 应 区 间 上 的 实 
数 根 。 

每 个 阵 元 只 接收 自身 发 射 的 信号 ,第 个 阵 元 
接收 的 目标 回 波 为 

R,(1)=AC—r)exp(j2rnf,t—r)), (6) 

其 中 ,z 为 双 程 时 延 ， 


(7) 


Ci C2 
将 接收 的 回 波 信号 与 exp(j2xf,t) 进 行 混 频 ,得 到 基 
带 信 和 号 ， 

B, (rt) =exp(—j2nxf,T)。 (8) 


2 点 扩展 函数 与 BP 成 像 算法 原理 


为 分 析 成 像 性 能 ,将 待 成 像 区 域 以 (Arz,Az ) 为 
单元 进行 网 格 划 分 ,通过 计算 网 格 点 与 各 阵 元 间 的 距 
离 ,构造 不 同 的 相位 补偿 。 成 像 区 域 网 格 划分 示意 图 
如 图 2 所 示 。 

将 目标 所 在 区 域 划分 为 NXN;, 个 网 格 点 , 任 
意 网 格 点 Q(x;,z; ) 对 第 n 个 阵 元 的 相位 补偿 项 为 


j2r 记 ( 双 二 a 


(9) 


Gi 二 exp{j2r 记 rr = exp 


图 2 网 格 划分 示意 图 


则 对 于 第 ”个 阵 元 ,整个 待 成 像 网 格 区 域 的 相位 补 
偿 和 矩阵 为 


91.1 1.2 OLN 
Sa 
ON, 外 ON a ON oN, 


将 所 得 待 成 像 区 域 相位 补偿 矩阵 对 点 目标 P 的 
基带 回 波 信号 进行 相位 补偿 ,可 得 第 n 个 阵 元 关于 
点 由 标 了 在 成 像 区 域 的 能 量 分 布 函 数 ; 

3 E,(P)=B,0,. (11) 
将 所 有 阵 元 在 该 区 域 的 能 量 进行 相干 累加 ,可 得 该 目 
标 虽 在 指定 区 域内 的 点 扩展 函数 ; 

~> Fs(P) = >E, CP). (12) 
对 沃 可 视 为 由 多 个 理想 散射 点 构成 的 目标 体 ,其 后 向 
投影 算法 得 到 的 目标 体 成 像 结 果 也 可 视 为 多 个 理想 
散 盘 点 各 自 的 点 扩展 函数 在 该 区 域内 的 春 加 : 


P 
OO I 一 >)F CD)。 (13) 
p=1 


3 ”介质 参数 估计 方法 


介质 的 参数 估计 是 地 下 目标 定位 和 成 像 的 关键 。 
只 有 对 介质 参数 做 出 较为 准确 的 估计 ,才能 使 偏 移 成 
像 获 得 较 好 的 结果 。 根 据 电 磁场 相关 理论 "” ,介质 
的 介 电 常数 es、 电导 率 c、 磁 导 率 ww 和 电磁 波 的 频率 w 
是 影响 电磁 波 传播 相 速 的 主要 因素 ， 


vy 二 L (14) 


| I+ 人 2) 得 
通常 情况 下 ,介质 多 为 磁 无 关 介 质 , 故 磁 导 率 mw 
可 认为 是 一 个 常数 。 分 析 式 (14) 可 知 , 电 导 率 为 0 
时 ,电磁 波 在 介质 中 的 传播 速度 与 频率 无 关 , 超 宽带 
信号 可 较为 清晰 地 进行 成 像 * 站 ;电导 率 不 为 0 时 ， 
不 同 频率 分 量 电磁 波 的 传播 相 速 不 同 , 产生 色 散 效 
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CilliidAiVO TFT 
应 ,进而 引起 波形 畸变 ,导致 测量 结果 出 现 误 差 。 而 
频率 分 集 阵列 各 阵 元 发 射 单 频 信号 ,可 对 不 同 频率 成 
分 单独 计算 其 相 速 。 

将 式 (14) 代 入 式 (9) ,可 得 : 


]2r 亡 Ee 1+( 了 ) + 下 
WnE1 


n 


2 
2R | 1+( 宇 ) | | | ， (15) 


当 电导 率 为 0 时 , 式 (15) 可 简化 为 
Di = exp(j2rf, (R33; Vuel + 2Rs,.i,; Vues ))。 
(16) 
从 式 (15)、(16) 可 看 出 ,介质 的 电 特 性 参数 \ 信 号 
的 频率 等 都 会 对 点 扩展 函数 的 成 像 质 量 产 生 一 定 影 
响 。 为 了 评价 成 像 的 质量 ,引入 图 像 箭 的 概念 。 
炉 是 某 一 体系 的 混乱 程度 的 表征 ,在 图 像 处 理 领 
域 , 图 像 箭 是 一 个 非常 重要 的 评价 指标 ,其 值 越 小 , 表 
明 图 像 的 聚焦 性 能 越 好 ,对 焦 越 准确 “"。 将 得 到 的 
各 离散 网 格 点 的 点 扩展 函数 灰 度 化 , 灰 度 值 为 工 , 则 
点 扩展 函数 灰 度 化 后 的 数值 可 用 如 下 矩阵 形式 表示 : 


11,1 11,2 ”ZN 


ij 一 EXp 


2 (17) 


IM,]1 M2 “”” MN 


设 不 同 灰 度 值 的 出 现 概率 为 p, , 则 图 像 粹 为 
on pi1 关 0。 (18) 


对 于 存在 目标 的 同一 块 待 成 像 区 域 , 当 采 用 不 同 的 介 
质 电 特性 参数 分 别 进行 成 像 时 ,图 像 的 聚焦 度 越 好 ， 
则 图 像 简 的 值 就 越 小 ,进而 对 于 介质 的 参数 ,存在 如 
下 估计 方法 : 

1) 通 过 工程 背景 等 先 验 信息 ,首先 确定 待 成 像 区 
域 介 质 的 相对 介 电 常数 和 电导 率 的 可 能 区 间 范 转 
[Le anyersnad,Lanunyano。 

2) 将 电导 率 分 别 设 置 为 其 可 能 区 间 范 围 的 两 端 
值 ,然后 设置 相对 介 电 常数 的 初始 值 ,并 选择 合适 的 
步 长 As,。 

3) 将 选 定 的 参数 值 进行 成 像 处 理 ,并 求 得 其 成 像 
结果 的 图 像 炉 。 

4) 改 变相 对 介 电 稼 数 的 值 ,s,=s, 十 As,, 重 复 步 
台 3) ,直至 相对 介 电 常数 变化 到 区 间 最 大 值 。 

5) 通 过 上 述 步骤 ,得 到 电导 率 分 别 取 端点 值 时 不 
同 s, 对 应 的 图 像 炉 值 ,对 其 取 平 均 ,图 像 炉 取 极 小 值 
时 所 对 应 的 s, 即 为 其 估计 值 。 

6) 选 择 合适 的 步 长 Ac ,将 初 值 设置 为 电导 率 的 
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估计 值 ,并 转 步 又 3)。 
至 电导 率 变化 至 区 间 最 大 值 。 


8) 同 理 , 图 像 炉 取 极 小 值 时 对 应 的 co 即 为 电导 率 
的 估计 值 。 


4 仿真 验证 及 结果 分 析 


4.1 频率 分 集 探 地 雷达 成 像 


选用 的 仿真 软件 为 MATLAB 2018a 和 
gprMax3D。 由 前 述 分 析 可 知 , 超 宽带 信号 在 介 
导 率 较 大 时 ,会 产生 色散 效应 , 且 会 影响 成 像 结果 
此 设计 了 以 下 实验 来 进行 验证 。 

实验 1 场景 如 图 3 所 示 , 场 景 大 小 为 0.6 mX 
036=m 的 电介质 , 蓝 色 实 线 为 阵列 所 在 位 置 ,在 场景 
中 窟 (0.3 m,0.3 m) 处 存在 一 金属 薄片 ,其 长 度 为 
2Cm ,用 以 模拟 理想 点 目标 。 超 宽带 探 地 雷达 及 频 
“名 和 集 蓉 地 埋 达 各 参数 如 表 1 所 示 。 


© 
CN > 介质 二 
©O 
CO 阵列 
GN 位 置 
© i 
ON - 目标 
之 
人 
SC 
有 图 3 实验 1 场景 示意 图 
© 
表 1 仿真 参数 
参数 妇 宽 带 GPR 频率 分 集 GPR 
中 心 频率 fo/GHz 1 1 
频率 增 量 Af/MHz 一 25.8 
阵 元 6 测 点 ) 个 数 N 31 31 
阵 元 ( 测 点 ) 间 隔 d/m 0. 02 0. 02 
阵列 ( 测 线 ) 长 度 工 /m 0.6 0.6 
介质 相对 介 电 常数 s, 6 6 
介质 电导 率 cl/(S,m 1) 0 0 
介质 电导 率 cv/(S,m 1) 0. 25 0. 25 
分 别 设置 介质 的 电导 率 为 0.0.25 S/m, 则 超 宽 


带 GPR 和 频率 分 集 GPR 对 不 同 电 导 率 的 介质 成 像 
结果 分 别 如 图 4、5 所 示 , 其 中 目标 真实 位 置 由 红色 圆 
图 标注 。 


为 便于 对 比分 析 , 将 成 像 结 果 均 进行 了 归 一 化 处 
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(b) 电导 率 为 0 时 频率 分 集 探 地 雷达 成 像 结果 
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(b) 电导 率 为 0.25 S/m 时 频率 分 集 探 地 雷达 成 像 结果 
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理 。 从 图 4、5 可 看 出 ,在 电导 率 为 0 时 ,2 种 模式 下 
的 探 地 雷达 均 能 够 较 好 地 对 目标 进行 定位 成 像 , 而 当 
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电导 率 设 置 为 0.25 S/m 时 , 因 超 宽带 信号 存在 色散 
效应 ,导致 其 成 像 结 果 存 在 一 定 误 差 , 而 随机 置换 频 
率 分 集 阵 列 , 仍 可 对 目标 进行 较为 精确 的 定位 成 像 。 


4.2 介质 参数 的 估计 


通过 上 述 实验 可 知 , 超 宽带 探 地 雷达 在 介质 存在 
一 定 电导 率 时 ,会 产生 信和 号 失真 ,使 得 恢复 的 目标 位 
置 不 够 精确 ,而 频率 分 集 阵 列 探 地 雷达 则 能 够 准确 地 
实现 目标 定位 。 当 目标 可 视 为 细 长 埋 地 目标 体 时 , 设 
计 以 下 实验 来 验证 本 方法 对 介质 参数 估计 值 的 准 
确 性 。 

实验 2 场景 如 图 6 所 示 , 设 置 为 1.2 mX1.3 m 
的 空间 范围 ,采用 阵 元 数 为 51 的 线性 阵列 ,其 距 下 层 
介质 0.3 m, 阵 元 间距 为 0.02 m, 初始 载 频 为 
TIGHz, 带 宽 为 200 MHz。 上 下 层 介 质 电 特 性 参数 如 


所 示 。 在 下 层 介质 中 心 位 置 (0. 6 m,0. 8 m) 处 ， 
一 细 长 理想 导体 圆柱 。 
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图 6 实验 2 场景 示意 图 


表 2 实验 2 各 介质 电 特 性 参数 


介质 ”相对 介 电 常 数 e。 电导 率 o/(S.m 1) ”相对 磁 导 率 y， 
空气 下 和 0 1.0 
粘土 6.0 0. 001 pe， 


采用 本 方法 ,根据 工程 背景 等 先 验 信息 ,确定 介 
质 的 相对 介 电 常数 及 电导 率 的 大 概 范围 ,并 选择 合适 
的 步 长 进行 多 次 迭代 ,得 到 相对 介 电 常 数 及 电导 率 的 
估计 值 ,如 图 7 所 示 。 

从 图 7 可 看 出 ,相对 介 电 常数 的 估计 值 为 6. 5， 
电导 率 的 估计 值 为 0.005 S/m。 采 用 预 设 的 实际 值 
和 得 到 的 介质 参数 佑 计 值 分 别 对 目标 进行 成 像 ,得 到 
的 成 像 结果 如 图 8 所 示 ,其 中 目标 真实 位 置 由 红色 加 
圈 标 注 。 


帅 等 :基于 频率 分 集 探 地 雷达 的 介质 参数 估计 方法 “，， 下 5 
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(a) 相对 介 电 常数 的 估计 结果 


(b) 已 导 率 的 估计 结 
图 7 对 介质 电 特性 参数 的 估计 


0.4 0.6 0.8 
方位 /m 
(a) 真实 参数 成 像 结 果 


0.4 0.6 0.8 


方位 /m 
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图 8 目标 成 像 结 果 图 


从 图 8 可 看 出 ,估计 值 和 真实 值 的 成 像 结果 存在 
一 定 误差 ,这 种 误差 主要 由 2 个 方面 引入 ,一 方面 ,由 
于 采用 电磁 仿真 软件 进行 网 格 离散 化 时 引入 了 一 定 
的 舍 入 误差 ,从 而 引起 电磁 波 传 播 速度 的 微小 变化 ; 
男 一 方面 ,在 计算 电磁 波 的 折射 点 位 置 时 也 会 引入 一 
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定 的 误差 ,致使 最 终 参 数 的 估计 值 及 相应 的 成 像 结 
虽 存 在 一 定 误 差 ,但 都 在 合理 范围 内 。 
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提出 了 一 种 基于 图 像 焙 的 地 下 介质 电 特 性 参数 
估计 方法 。 推 导 了 频率 分 集 阵列 的 点 扩展 函数 ,并 对 
比 了 点 扩展 函数 与 后 向 投影 算法 之 间 的 异同 ;对 超 宽 
人 带 信 号 在 有 耗 介 质 中 的 色散 现象 进行 了 分 析 ,对 不 同 
介质 条 件 下 的 超 宽 带 探 地 雷达 与 频率 分 集 探 地 雷达 
的 成 像 结 果 进 行 了 对 比 。 实 验 结果 表明 ,频率 分 集 探 
地 雷达 比 传统 超 宽 带 探 地 雷达 在 目标 定位 成 像 方面 
具有 更 好 的 效果 ,在 介质 有 耗 的 情况 下 ,依然 能 够 较 
为 精确 地 对 目标 进行 定位 成 像 。 因 此 ,本 方法 能 够 在 
合理 的 误差 范围 内 ,对 地 下 介质 参数 进行 有 效 估计 。 
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